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We present a new approach to quantifying pole parameters of single-channel processes based on 
Laurent expansion of partial wave T-matrices. Instead of guessing the analytical form of non-singular 
part of Laurent expansion as it is usually done, we represent it by the convergent series of Pietarinen 
functions. As the analytic structure of non-singular term is usually very well known (physical cuts 
with branhcpoints at inelastic thresholds, and unphysical cuts in the negative energy plane), we show 
that we need one Pietarinen series per cut, and the number of terms in each Pietarinen series is 
automatically determined by the quality of the fit. The method is tested on a toy model constructed 
from two known poles, various background terms, and two physical cuts, and shown to be robust 
and confident up to three Pietarinen series. We also apply this method to Zagreb CMB amplitudes 
for the N(1535) 1/2- resonance, and confirm the full success of the method on realistic data. This 
formalism can also be used for fitting experimental data, and the procedure is very similar as when 
Breit-Wigner functions are used, but with one modification: Laurent expansion with Pietarinen 
series is replacing the standard Breit-Wigner T-matrix form. 

PACS numbers: 11.55.-m, 11.55.Fv, 14.20.Gk, 25.40.Ny. 
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The recent Camogli workshop [Ij has finally inaugu- 
rated the fact that poles, and not Breit-Wigner param- 
eters determine and quantify resonance properties, and 
that they should be used a link between scattering the- 
ory and QCD. However, at the same time, the ques- 
tion of finding an adequate procedure to extract them 
from single-channel T-matrices has been opened. Exper- 
imentalists are quite familiar with fitting the data with 
Breit-Wigner functions (either with constant parameters 
and very general backgrounds, or with energy dependent 
mass and width), but are inexperienced when actual com- 
plex energy poles have to be used. A simple procedure 
for pole extraction is still missing. Up to now poles have 
usually been extracted from theoretical single or multi- 
channel models fitted to the data using either of stan- 
dard pole extraction methods: analytic continuation of 
the model functions into the complex energy plane [JHS] , 
speed plot ^, time delay [8j, N/D method [5], regular- 
ization procedure [TD], etc. However, this required solv- 
ing standardly elaborated single/coupled-channel models 
and analyzing the obtained analytic solution which im- 
plicitly contained both parts: singular and background. 
Hence, the analytic form of the full solution varied from 
model to model; pole-background separation was not 
unique, and this introduced uncertainties in pole extrac- 
tion procedures. The intention of this paper is to offer 
simple, robust and confident method how to obtain scat- 
tering matrix poles, but maximally avoid the need to 



speculate about type and form of background terms. We 
base our analysis on Laurent expansion of partial wave 
T-matrices which uniquely separates singular from finite 
terms, and treat singular and finite terms separately. Our 
main assumption is that all scattering matrix poles are 
of the first order. 

The starting point of our method is Laurent expansion 
for a function with one pole: 



T{u) 



+ B^iLo); a_i,a;o,a; e C. (1) 



where a_i and luq are residuum and pole position re- 
spectively, and function B^{uS) = Yl^=o°'n(!^o ~ ^)'^ is 
regular everywhere in the complex energy plane. How- 
ever, functions we analyze in reality may and do contain 
other poles for w 7^ wq, so if we iterate this procedure we 
can without loss of generality write down: 

k {i) 

T{uj) = V— ^+B^(cj); J'\,oj,,uj eC. (2) 

— ^ L0{ — iJ 

2—1 

where k is number of poles. a_[ and Wi are residua and 
pole positions for i-th pole respectively, and B^{uj) is a 
function regular in all uj ^ uji. 

This approach has been already investigated, but the 
freedom in choosing the exact analytic form of the back- 
ground contribution B^{uj) has been introducing severe 
ambiguities. 
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TABLE I: Toy model parameters and fitted parameters. Input parameters are given in boldface, and results of a fit in normal 
font. Table is given in GeV units, and Fj = —2 Wi. 



Ci 






B2 


n gi Ml Fi r2 g2 M2 F2 


a xp Ni 


13 XQ N2 


7 xr N3 




Toy-model 






0.1 0.09 1.65 0.165 0.09 0.06 2.25 0.2 










Fitted results 




Strategy a. 




1 











0.100 0.089 1.649 0.165 0.090 0.060 2.249 0.200 


2.48 0.97 5 






0.03 





1 








0.099 0.090 1.650 0.165 0.090 0.060 2.249 0.199 


3.97 3.97 5 






0.01 








1 


1 


0.098 0.091 1.650 0.165 0.090 0.060 2.250 0.200 


1.19 -14.94 7 






0.2 








-1 


-1 


0.099 0.089 1.649 0.1649 0.089 0.059 2.249 0.199 


0.99 -9.63 7 






0.01 


1 





1 


1 


0.103 0.100 1.653 0.171 0.101 0.067 2.249 0.221 


0.71 -0.23 11 






28 


1 





1 


1 


0.099 0.090 1.650 0.164 0.089 0.060 2.250 0.199 


-2.04 -17.58 5 


4.27 0.97 5 




0.28 


1 





- 1 


- 1 


0.097 0.087 1.651 0.161 0.090 0.060 2.250 0.201 


0.90 -0.39 20 






22.0 


1 





-1 


-1 


0.099 0.089 1.649 0.164 0.090 0.059 2.249 0.199 


2.96 -8.97 6 


1.56 0.97 6 




1.00 





1 


1 


1 


0.107 0.088 1.646 0.166 0.093 0.048 2.239 0.197 


2.06 -0.89 10 






114.79 





1 


1 


1 


0.099 0.090 1.650 0.165 0.090 0.060 2.250 0.200 


1.94 -16.33 5 


6.42 3.97 5 




0.02 





1 


-1 


-1 


0.090 0.086 1.651 0.156 0.095 0.058 2.248 0.202 


0.969 -0.37 12 






238.38 





1 


-1 


-1 


0.099 0.090 1.650 0.165 0.090 0.060 2.250 0.200 


0.81 -7.89 8 


1.24 3.97 8 




0.06 


1 


1 


1 


1 


0.085 0.102 1.663 0.171 0.087 0.075 2.262 0.216 


1.09 -2.64 10 






328.19 


1 


1 


1 


1 


0.098 0.086 1.650 0.161 0.095 0.058 2.247 0.199 


0.44 -0.47 9 


1.95 3.97 8 




70.37 


1 


1 


1 


1 


0.099 0.090 1.650 0.164 0.089 0.061 2.251 0.200 


4.19 -22.99 5 


2.22 3.98 5 


1.67 0.97 3 


0.24 


1 


1 


-1 


-1 


0.090 0.105 1.657 0.182 0.078 0.061 2.260 0.189 


1.38 -3.12 10 






467.54 


1 


1 


-1 


-1 


0.095 0.098 1.654 0.173 0.086 0.061 2.254 0.198 


0.61 -0.20 9 


25.91 3.98 8 




60.94 


1 


1 


-1 


-1 


0.100 0.090 1.650 0.165 0.090 0.060 2.250 0.200 


1.85 -6.25 3 


16.36 3.97 3 


1.32 0.98 3 


0.72 


Strategy b. 




1 











0.069 -0.111 1.647 0.155 0.081 0.055 2.252 0.193 


0.96 0.994 18 






0.83 





1 








-0.101 -0.084 1.649 0.165 -0.088 -0.061 2.249 0.200 


0.90 3.94 9 






0.33 








1 


1 


0.329 0.247 1.649 0.165 -0.077 0.076 2.249 0.200 


1.00 -0.76 8 






0.01 








-1 


- 1 


0.160 0.158 1.649 0.165 0.235 0.109 2.249 0.200 


0.65 -8.41 8 






0.01 


1 





1 


1 


0.114 -0.049 1.657 0.156 -0.146 0.418 2.254 0.187 


0.52 0.00 11 






12.3 


1 





1 


1 


-0.452 -0.036 1.649 0.165 0.096 -0.064 2.249 0.200 


1.00 -5.81 5 


2.06 0.97 5 




0.01 


1 





-1 


-1 


-0.226 0.160 1.645 0.166 0.017 0.110 2.249 0.208 


0.34 -0.02 20 






7.12 


1 





-1 


-1 


0.116 -0.010 1.650 0.164 0.226 -0.202 2.249 0.199 


2.14 -0.18 6 


1.70 0.98 6 




0.03 





1 


1 


1 


0.320 0.060 1.643 0.166 0.036 0.111 2.244 0.229 


1.00 -0.31 10 






43.59 





1 


1 


1 


-0.096 -0.092 1.650 0.165 -0.090 -0.059 2.250 0.200 


1.11 -3.81 5 


1.60 3.97 5 




0.01 





1 


-1 


-1 


0.062 0.143 1.653 0.184 0.202 0.329 2.268 0.225 


0.85 -0.05 9 






102.28 





1 


-1 


-1 


0.101 0.291 1.650 0.165 0.090 0.062 2.250 0.199 


4.25 -66.59 6 


20.03 3.97 6 




0.01 


1 


1 


1 


1 


0.052 -0.092 1.656 0.143 0.239 0.282 2.235 0.179 


1.18 -0.25 12 






46.02 


1 


1 


1 


1 


-0.058 -0.122 1.662 0.167 -0.054 -0.081 2.258 0.185 


0.71 -0.72 6 


1.38 4.00 6 




25.22 


1 


1 


1 


1 


-0.318 0.258 1.648 0.165 0.073 -0.104 2.247 0.207 


0.67 -9.20 7 


0.17 3.99 7 


0.18 1.00 7 


1.09 


1 


1 


-1 


-1 


0.098 0.365 1.640 0.164 0.126 0.061 2.247 0.118 


1.66 -1.16 12 






25.63 


1 


1 


-1 


-1 


-0.012 0.375 1.649 0.165 0.080 0.077 2.251 0.200 


1.36 0.92 8 


2.62 3.98 7 




0.53 


1 


1 


-1 


-1 


-0.088 0.069 1.65 0.164 -0.015 0.370 2.249 0.201 


1.11 -1.62 7 


3.53 4.02 7 


1.56 0.97 5 


0.23 
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The novelty of our approach is that we propose to 
avoid discussing the arbitrariness in aU possible choices 
for the background function B^{uj) by replacing it 
with Pietarinen expansions in power series using a com- 
plete set of functions with well known analytic properties. 

Pietarinen series 

If F{uj) is a general, unknown analytic function having 
a cut starting at co — xp, then it can be represented in 
a power series of Pietarinen functions in the following 
way: 

N 

F{uj) = ^c„Z(w)", weC 

ri=0 

a — Jxp — Lo 
Z(u}) = — — , Cn,xp,a eK, (3) 

a + ^xp - uj 

with the a and c„ being tuning parameter and coeffi- 
cients of Pietarinen function Z(uj) respectively. 

The Pietarinen series have been proposed and in- 
troduced by Ciulli [11] and Pietarinen [12], and have 
been with great success used in Karlsruhe-Helsinki 
partial wave analysis [13] when invariant amplitudes 
have been expanded in as many as 50 terms. The essence 
of the approach is the fact that a set n = 1, oo) 

forms a complete set of functions defined on the unit 
circle in the complex energy plane having branch cut 
starting at w = xp, but the analytic form of the function 
is at the beginning yet undefined. The final form of the 
analytic function F{u;) is obtained by introducing the 
rapidly convergent power series with real coefficients, 
and the degree of the expansion is automatically deter- 
mined in fitting the input data. 

The application of Pietarinen series to scattering theory 

The analytic structure of each partial wave is a well 
known fact. Each partial wave contains poles which 
parameterize resonant contributions, it has cuts in the 
physical region starting at thresholds of elastic and 
all possible inelastic channels, and finally there are 
t-channel, u-channel and nucleon exchange contributions 
quantified with corresponding negative energy cuts. 
However, explicit analytic form of each cut contribution 
is not known. Instead of guessing the exact analytic 
form of all of them, we propose to use one Pietarinen 
series to represent each cut, and the number of terms 
in Pietarinen series will be determined by the quality 
of fit to the input data. So, in principle we have one 
Pietarinen series per cut, the branch-points xp,xq... are 
known from physics, and coefficients are determined by 
fitting the input data coming from real physical process. 
In practice, we have too many cuts (especially in the 
negative energy range), so we reduce their number by 
dividing them in two categories: all negative energy 
cuts are approximated with only one, effective negative 
energy cut represented with one Pietarinen (we usually 
denote its branchpoint as ccp), and each physical cut 



is represented with its own Pietarinen series with 
branch-points determined by the physics of the process 

{xq,Xr...). 

So, the set of equations which define Laurent expan- 
sion-|-Pietarinen series method (L-l-P method) is: 

k (i) 

M N 

B'^iuj) = ^ c„ Z(cj)" + ^ d„ T^(w)" + • • • 

n=0 n=0 

a + ^/xp - uo fj + ^XQ~u! 

a^i, Wi, a; G C 

Cn,xp,dn,XQ,a, (3... e M 

and k,M,N...e N. (4) 

As our input data are on the real axes, the fit is per- 
formed only on this dense subset of the complex energy 
plane. All Pietarinen parameters in set of equations Q 
are determined by the fit. 

Let us observe that the class of input functions which 
are convenient to be analyzed with this method is 
quite wide. One may either fit partial wave amplitudes 
obtained from any of theoretical models, or even ex- 
perimental data directly. In either case, T-matrix is 
represented by the set of equations Q, minimization 
function is defined (usually type), and fitting is done. 

Let us summarize our fitting procedure 

First of all let us observe that in the strict spirit of 
the method, physical Pietarinen branch points xq, xp, ... 
should not be fitting parameters. Namely, as we have 
declared that each known cut should be represented by 
its own Pietarinen series having the analytic structure 
of the very cut, Pietarinen branch points should be 
fixed to known physical branch points. And this indeed 
is so in the ideal case. However, in realistic case the 
situation is very similar to the background situation. 
Namely, we can never include all physical cuts in the 
multi-channel process, so instead of taking them all, we 
represent them by a smaller subset. So, for physical 
energy range too, Pietarinen branch points XQ,xp,... 
are not constants; we have to relax them and allow them 
to vary as fitting parameters. Later on in this paper we 
shall demonstrate that when we do it, physical branch 
points still naturally converge towards branch-points 
which belong to channels which dominate certain partial 
wave, but do not actually correspond to them exactly. 
The proximity of the fit results to exact physical branch 
points describes the "goodness of the fit, namely it 
tells us how well certain combination of thresholds 
indeed is approximating certain partial wave. And this, 
together with the choice of the degree of Pietarinen 
polynomial represents the model dependence of our 
method. We have, of course, never claimed that our 
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method is model independent because there is no such 
thing as model independence. However, the method 
fixes its model dependence naturally, by fitting to the 
real data. It chooses the simplest function with the given 
analytic properties which fit the data, and increases 
the complexity of the function only when the actual 
data require so. We shall later on demonstrate that for 
1/2- partial wave (Su) the fit chooses ttN elastic and 
rj production branch point, while for 1/2+ partial wave 
(Pii) instead of choosing r] production branch point, it 
settles close to the ttttN threshold which is dominant 
for this partial wave. Let us observe that we are still 
limited to the stable two-body channel representation 
(real branch point), so quasi stable two body channels 
like pN and aN are still not included. L-f P method can 
do it by allowing the branch points xq, xh, ... to become 
complex numbers, but it simply has not been done yet. 

We first start with minimal number of poles, the 
minimal number of Pietarinen functions (we choose 
only dominant inelastic channels), minimal number of 
fitting parameters a, /3,c„, dn and Pietarinen branch 
point xp,XQ,... being close to actual physical branch 
points. We usually start with N — 3. The reduced 
is analyzed, and the quality of the fit is visually 
inspected by comparing fitting function with fitted data. 
If the fit is unsatisfactory (reduced is high, or fit 
visually does not reproduce fitted data), the number of 
Pietarinen parameters c„, d„ ... is increased by one. The 
fit is repeated, and the quality of the fit is re-estimated. 
This procedure is continued until we have reached the 
sufficient number of Pietarinen terms so that we are able 
to reproduce the fitted data. If the quality of the fit is 
still unsatisfactory, we first increase the number of poles 
and repeat the procedure. If no improvement is achieved 
by increasing the number of poles, then we increase 
the number of Pietarinen functions by one taking the 
next branch-point (next physical threshold), and re- 
peat the procedure until the agreement has been reached. 

Testing the method 

We have tested the method in three ways: 

I) We have constructed a toy-model function imi- 
tating the physical reality as close as possible (known 
pole parameters, two positive energy cuts and various 
background contributions), and verified how well our 
method reproduces the input parameters; 

and to see how our method works in realistic cases 

II) We have used our L-t-P method to extract pole 
parameters from single-channel amplitudes of a known, 
published model (N(1535) 1/2- amplitude of Zagreb 
CMB model [l^) and compared the L-l-P result with 
pole parameters from the original publication, and 

III) We have used our L+P method for something 



what could never be done before, we have extract pole 
positions from partial wave data of ref. |15j without 
assuming any final functional form for the scattering 
matrix regular part. 

I) Fitting the toy model 

In principle, we should have defined a toy model 
input data T*y{ijjj) by defining a toy model function and 
by normally distributing its values in order to simulate 
the statistical nature of real measured data. However, as 
the main goal of this paper is to establish the validity of 
the approach, we have restricted our analysis to infinitely 
precise data by using non-distributed toy-function values, 
and using the statistical weight Wj of 5 %. This enables 
us to test the details of our theoretical assumptions, but 
gives unrealistically low values. Testing the capac- 
ity and limitations of the L-fP method with realistic 
data for the toy model is deferred to another publication. 

A toy model function is constructed by assuming 
a typical analytic structure of a partial wave: it is 
constructed as a sum of two poles, two physical cuts and 
several non-resonant background contributions. The 
function representing physical cuts is constructed from a 
function 5fte(x,a) = yjx'^ — Aax/2x having a cut starting 
from xq = 4a on the real axes^, and the analyticity is 
imposed through the once subtracted dispersion relation 
$(x, a) = ^ f °° , dx'. However, to sim- 

plify the demonstration of usability of L-l-P method, we 
replace all negative energy cuts with two poles deep in 
unphysical region. In spite of looking rather restrictive, 
such an approximation is fairly justified. Namely, we 
know that each cut can be represented by the infinite 
sum of poles, and as negative cut is indeed very far from 
the region of interest, replacing it with only two out of 
infinite number of poles is a good approximation (see 
Cutkosky CMB approach [4]). 
So our toy -model function is given as: 



T'y{u:) = 



M. 



ty , ty 

ty 



ty 



(5) 



$(iu, a) 
where 



Ci $(a;, 0.25) + C2 $(w, 1.) + B*^(w), 

^yuJ^^^4n^+UJ) 2a — uj — ^Loi—Aa + cj) 
m 



27rw 



2a 



10. 



-10. -uj-ib. 



Bo 



10. 



-6.-a;-z4.' 



rf,gf,Mf,W,*^e 



^ The type of the function used for physical cuts comes from the 
phase space factor for two body reactions 4>{s) = \J A(s) /2s, with 
A(s) = - 2s{M^ + ml) + (M^ - m^)^, and taking m,r = M. 
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<n[GeV] (l)[GeV] 
(a) (b) 



FIG. 1: Toy model function. In figures a) and b) we give 
toy-model data for the function with two poles and two cuts 
for two choices of background parameters, for {Bi — 1, B2 
= 1) and {Bi = -1, B2 = -1) respectively. Dashed-dotted, 
dashed and full lines PI, P2 and P3 give the quality of the fit 
for solutions with one, two and three Pietarinen expansions 
respectively. 

Toy model parameters for all test cases are chosen to 
resemble physical reality as much as possible, and are 
given in Table |l] with bold face characters. 

We have applied two fitting strategies: 

a) fitting both, real and imaginary part of the toy- 
model data (imitating the physical situation when a 
complete experiment is performed, and the full T-matrix 
is known); and 

b) fitting only absolute value of the toy-model function 
(imitating physical situation when only incomplete data 
set is available, and these are usually differential cross 
section data). 

Minimization function for case a): 
= Xf + AxPe„; 

Xf = 5] |r*n^,)-n^.)lVwf. (6) 

Minimization function for case b): 

= Xf + ^XPen, 

Xf = 5]||T*na;,)p-|T(c.,)n/^'- (7) 

In both cases Wj is corresponding statistical weight and 
XPen = X^feLi is Pietarinen series penalty function 
|12j which guarantees that the number of power-series 



terms is minimal. Coefficient A is determined in a fit, and 
serves to match the size of Xf versus XPen contributions. 

In situation a) we expect a unique solution, and in sit- 
uation b) we expect a full set of solutions because the 
relative phases of the fitted functions are still undeter- 
mined. 

Let us introduce a notation [a,b,c,d] meaning by defi- 
nition: (Ci= a, 6*2= b, Bi— c, i?2= d). 

We have tested the validity of the model for two back- 
grounds labeled A=:[a,b,l,l] and B=[a,b,-1,-1] which 
strongly contribute to the toy-function strength in the 
observed resonance range (1 < w < 3 GeV), and produce 
drastically different form of the toy-function absolute 
value [see Fig. ([T])]. We have decided to characterize the 
type of the background according to the absolute value 
shape (what basically corresponds to the differential 
cross section): while background A produces typical 
"two peak" resonance structure, background B produces 
a very atypical behavior in the first-second resonance 
region. If these numbers were differential cross sections 
extracted from experiment, one could not easily say 
whether they in case of background B represent a 
resonance, or some other interference effect. Therefore, 
we believe that these two backgrounds are the worst 
case scenario for the L-|-P method to separate resonance 
and background contributions. 

Results for fitting strategy a) 



With toy model data we tested the functionality 
of the model for various combinations of toy-model 
ingredients. The aim is to verify basic concepts of 
L+P method for the ideal data set. All results of the 
toy-model L-l-P fit are shown in Table |T] Toy model 
function with all ingredients included, and for two 
different background contributions {Bi = 1, B2 = 1) 
and {Bi = -1, B2 = -1) is depicted in Figs, p] a. and 
b.). 

Two poles, no background one cut 

We fit the toy model data generated by the toy function 
with [1,0,0,0] and [0,1,0,0]. We establish that one 
Pietarinen expansion in L-|-P formalism is sufficient to 
reproduce the input data (see low value). However, 
Pietarinen branch point xp differs for both solutions (1 
and 4) indicating that the cut structure is reproduced. 
Both, residua and pole positions are perfectly repro- 
duced. 

Two poles, background, no cuts 

We fit the toy model data generated by the toy function 
with [0,0,1,1] and [0,0,-1,-1]. We establish that one 
Pietarinen expansion in L+P formalism with cut at 
Xp ^ is sufficient to reproduce the input data (see 
low x^ value). Pietarinen branch point xp differs for 
both solutions indicating that the different backgrounds 
are fitted. Both, residua and pole positions are perfectly 
reproduced. 

Two poles first cut, first background 

We fit the toy model data generated by the toy function 
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with [1,0,1,1] (background A). We see that for one 
Pietarinen expansion reduced value is high, and the 
Pietarinen branch point is negative. To improve the fit 
we had to introduce second Pietarinen series. The fit is 
improved, and the fitting has converged. So, in spite of 
the fact that we have represented the negative cut as a 
sum of unphysical poles, we still see it as a cut, so we 
need two Pietarinen expansions to take both cuts into 
account. While in case of one Pietarinen the branchpoint 
is negative, in case of two Pietarinens one branch point 
is negative, and second one converges towards 1 what is 
the toy- function branch point. Both, residua and pole 
positions are perfectly reproduced. 
Two poles first cut, second background 
We fit the toy model data generated by the toy function 
with [1,0,-1,-1] (background B). The situation of as 
for first background [1,0,1,1] is reproduced. We again 
need two Pietarinen series. Overall conclusion is that 
regardless of the type of background, L-fP formalism 
works. Both, residua and pole positions are perfectly 
reproduced. 

Two poles second cut, both backgrounds 
We fit the toy model data generated by the toy function 
with [0,1,1,1] and [0,1,-1,-1]. Situation is exactly repro- 
duced as for first cut, but with one difference: Pietarinen 
branchpoint xp for two Pietarinen expansions converges 
towards xp = 4 what is exactly second toy- model 
cut. Both, residua and pole positions are perfectly 
reproduced. 

Two poles two cuts, first background 
We fit the toy model data generated by the toy func- 
tion[l, 1,1,1] first with only one Pietarinen expansion, 
and afterwards with two. In both cases reduced 
value is, as expected, high. The value of is better 
for two Pietarinens, but the result is still unsatisfactory. 
Adding third Pietarinen fixes the problem. The fit 
has converged, and Pietarinen cuts are consistent with 
expectations: first one is negative, and second and 
third are close to 1 and 4, namely close to toy-model 
branch-points. Both, residua and pole positions are 
perfectly reproduced. 

Two poles two cuts, second background 
Conclusions are identical as for the first background. 
Consequently, the L-l-P formalism is invariant with re- 
spect to the relative size and phase between background 
and physical contributions. 

Results for fitting strategy b) 

All conclusions as for fitting strategy a) are repro- 
duced with one major exception: 

we only reproduce pole positions, residua are quite 
arbitrary. As a matter of fact, when fitting strategy b) 
was used, we have been able to produce a whole class 
of solutions with almost identical ^'^d very different 
residua, exactly what was to be expected. 







1 




0.4 




Data set 








Fit p3(2 poles) 








Fit p3(3 poles) 




0.2 

0.2 








1,1,1, 



W[MeV] 



' it ' ' 


^ ' , — 1,1 


Data set 

Fit p3(2 poles) 

Fit p3(3 poles) 

1 





oi r- \ < \ < 1— 

1000 1500 2000 2500 



W[MeV] 

FIG. 2: Zagreb Sll. Dashed and full lines gives the quality 
of the fit for two and three pole solutions in both cases with 
three Pietarinen series. 



//; Fitting ttN elastic N(1535) 1/2- 

To illustrate how usable L-l-P method is in reality, we 
have decided to fit one realistic set of input data. We 
have taken vrA^ elastic N(1535) 1/2- partial wave from 
Zagreb CMB model, and fitted it with L+P method. We 
have chosen fitting strategy a) and fitted both, real and 
imaginary part. We have to mention that we did not 
have much use of the Penalty function in fitting the toy- 
model, but for realistic data from Zagreb CMB using the 
penalty function was unavoidable. Results are shown in 
Table [n] and Fig. 

We know that ttTV elastic scattering has at least two 
branch-points in the physical region [t:N elastic thresh- 
old in Zagreb model at w = 1.076 GeV and 77 production 
In Zagreb model at w = 1.486 GeV), background con- 
tribution, and at least two poles. So we know that our 
L-l-P solution should have minimally two poles and three 
Pietarinen series. Therefore, we have started the fitting 
procedure with two poles, and three Pietarinen series, 
and obtained a reasonable solution. The solution was 
acceptable visually, and by x^- Two Pietarinen branch- 
points came out very close to physical thresholds, the 
third Pietarinen branch-point came out far in the neg- 
ative energy plane describing the background contribu- 
tion from the negative cut, and pole parameters have 
been quite close to parameters obtained in analytic con- 
tinuation. We wondered if the fit could be improved by 
increasing the number of poles by one, so we have tried 
a three pole-three Pietarinen fit. The reduced was 
somewhat improved, the pole parameters came some- 
what closer to the analytic continuation value, but in gen- 
eral the third resonance was uncertain and rather wide. 
It is interesting to observe that both physical Pietarinen 
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TABLE II: Parameters of the fit for the tvN elastic N(1535) 1/2- partial wave from Zagreb CMB model. Table is given in MeV 
units, and Fi = — 2 Wi. 





Ml Ti 


Ma Ta 


Ms Ta 


xp, Ni 


XQ, Na 


Xr, Na 




Zagreb CMB 


1521(14) 190(28) 


1646(8) 204(17) 


1790(26) 420(45) 






p... 2 poles 

SPictarinens 


1525 120 


1653 189 




-696, 15 


1058, 15 


1484, 15 


2.53 


p... 3 poles 

^ SPiotarinens 


1529 146 


1647 192 


1801 2321 


933, 15 


1057, 15 


1482, 15 


1.31 



branch-points remain where they were (close to physi- 
cal thresholds), but the third came out much higher. So, 
L-l-P method conclusively gives two resonances in Zagreb 
CMB amplitudes, but indicates the existence of the third 
one which is not well defined. 

We ended up with using three Pietarinen series 
with one threshold in unphysical range representing 
the Pietarinen series for background contributions, and 
two thresholds in the physical range. One was at 
xq — 1.057 Gev (near physical threshold), and sec- 
ond at Xr — 1.482 Gev (near rj production threshold). 
Pole positions of Zagreb model are almost perfectly 
reproduced. Using L-t-P method we obtained three 
poles: Ml = 1.529 - i 0.073, M2 = 1.647 - i 0.096 and 
= 1.801 — i 1.160 GeV, what is to be compared to Za- 



greb CMB poles from ref. [14]: M 
Mf» = 1.646 



Zg 



i 0.102 and M^<^ = 1.790 



1.521 -i 0.095, 
-i 0.210 GeV. 
First two poles agree almost ideally, the existence of third 
pole around 1.8 Gev is allowed, but single-channel data 
are insufficient to pin it down more precisely as expected 
from coupled-channel calculation. 

As Zagreb CMB analytic continuation was performed 
in three channel model, we conclude that two resonances 
is the best what a single- channel method like L-l-P can 
give. We also allow for some deviations in quantitative 
values of two well established resonances, since L-l-P is, 
again, only single-channel model. A correct recipe would 
be to repeat the L-f-P procedure on amplitudes from all 
channels, and then make an analysis. Extending the anal- 
ysis to other channels might shift pole position somewhat 
deeper in the complex energy plane, single channel anal- 
ysis as we did it now is, however, "satisfied" with poles 
being somewhat closer to the real axes. 

Ill) Poles from experiment: 
Fitting GWU-SES partial wave data 

The novelty of L+P method is that it allows extracting 
poles from data coming directly from experiment, i.e. to 
analyze the numbers which are obtained using only sta- 
tistical methods with very little underlying theory. We 
shall illustrate this feature on GWU single energy (SES), 
and energy dependent (SP06) solutions [15,. 

We have extracted pole positions from both, GWU- 
SES (" experiment" ) and energy dependent SP06 solution 
which is obtained from GWU-SES with the theoretical 
analysis based on using the polynomial form of Chew- 
Mandelstam K-matrix model, and compared results. 




Pll 




1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 2100 
(i)[McV] 



FIG. 3: GWU SP06 and SES. 



Let us point out that GWU group is able to extract 
poles from their SES only through creating their energy 
dependent SP06 solution, and analytically continuing it 
into complex energy plane. Our method allows us to fit 
both, SES and SP06, and compare results. 
We came to a very interesting conclusion: 
We can fit both, GWU-SES and SP06, by using two poles 
and three Pietarinen series, first threshold being in un- 
physical range, second near physical threshold and third 
near ttttN threshold (observe that Sn was fitted with 
Pietarinen near rj production threshold). Results of a 
fit are given in Table |III[ In boldface we give pole posi- 
tions from original publication. The original publication 
finds Roper poles on two Riemann sheets (on P[221] and 
on P[121]), and we give them both. However, our L-fP 
method is restricted to the first Riemann sheet only. 

From Table IIIII we see that for SP06 the second 
pole in our fit is completely undetermined, while for 
GWU-SES it is definitely established, and close to 1.7 
GeV. And that is very similar to what GWU group 
claims: they claim that they do not need Pii(1710) in 
their analysis at all. However, their analysis is based 
solely on energy dependent analysis. We agree with 
them in a sense that we also do not need Pii(1710) 
to fit their energy dependent solution SP06. How- 
ever, when fitting GWU-SES what they could not do, 
our fit REQUIRES the presence of the Pii(1710) 
state! Our results are: Aff™^ = 1.358 - i 0.094, 
A/|P*'6 = 2.694 - i 0.893 GeV for SP06 (sec- 
ond resonance completely undetermined), and 
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TABLE III: GWU SES and SP06 from ref. (15| . Poles from original publication are given in boldface. Table is given in GeV 
units, and Fi = — 2 Wi. 





n gi Ml Fi r2 32 M2 F2 


Q Xp Ni 


/3 XQ N2 


7 xr N3 


Xr 


GWU SES 






GWU SP06 


1.388 164 

1.358 162 - - - - 




L+P SES 


0.008 -0.059 1.358 0.156 -0.001 0.006 1.679 0.106 


1.109 -2.409 15 


0.573 1.074 15 


0.677 1.243 15 


1.815 


L+P SP06 


0.001 -0.065 1.357 0.188 0.180 0.504 2.694 1.786 


0.987 -1.195 15 


1.125 1.024 15 


0.711 1.208 15 


0.005 



Mps ^ 1^ 358 „ j 0.078, Mf^s ^ j^ gjg _ j o.053 GeV 
for SES (second resonance firmly established). 

So, our L-l-P method allows us to explicitly see that 
their SES contain Pii(1710), and their energy dependent 
analysis smeared it out. Results of a fit are given in 
Fig-i 

Let us also observe that taking only two Pietarinen 
series for physical range is still an approximation. 
Namely, in both fits given in Table |III| second and third 
threshold are close, but not exactly on top of first two 
physical thresholds. The first threshold for SES is very 
close to irN channel opening (at 1.076 GeV), while the 
second one is notably above ttttN threshold (at 1.218 
GeV) . For SP06 fit thresholds show opposite trend - first 
threshold is significantly lower, and second threshold 
is very close to ttttN value. That indicates the level 
of precision of our pole results, and stimulates further 
improvement of the model. However, even this level 
of precision is sufficient for pursuing the main argument . 

Conclusion 

We propose a new single-channel pole extraction 
method based on Laurent decomposition and Pietarinen 
series. Instead of guessing the unknown T-matrix 
functional form as it has been done up to now, we 
propose to use the Laurent's expansion of single-channel 
partial wave T-matrices to obtain their pole positions. 
We make use of a fairly well known analytic structure of 
each partial wave, and expand the unknown regular part 
of the Laurent expansion in power series of Pietarinen 
functions using one Pietarinen series for each known 
cut. Finally, we fit the input data and get the simplest 
function for the regular Laurent part which has the 
correct analytic structure and together with the known 
number of poles reproduces the input. The resulting 
pole parameters are the genuine poles of our unknown 
input function. Even to the surprise of the authors, the 
fitting procedure in analyzed test cases turns out to be 
extremely stable and reliable regardless of the number 
of Pietarinen series, and extracted pole parameters turn 
out to be fairly confident. We warn the reader that this 
is a single- channel method^ and it will recognize only 
those resonances which strongly couple to the observed 
channel. Other resonances will be only hinted at, and 
have to be pinned down by analyzing other channels. 
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Appendix 

Pietarinen expansion coefficients of solutions presented in Tabie^ (Strategy a). 



Startcgy 


a: 










One Pietarinen expansion 










Ci 


C2 


Bi 


B2 


CO 
Cll 


Cl 
Cl2 


C2 
Cl3 


C3 
Cl4 


C4 
Cl5 


C5 
Cl6 


C6 
Cl7 


C7 
Cl8 


C8 
Cl9 


eg 

C20 


eio 
621 


1 











1.079 


-1.080 


0.155 


-0.175 


0.044 


-0.022 


- 













1 








0.947 


-1.031 


0.205 


-0.153 


0.044 


-0.012 


- 


- 


- 


- 


- 








1 


1 


26.944 


62.508 


20.310 


-51.562 


-20.914 


55.819 


58.087 


17.230 


- 


- 


- 








-1 


-1 


-54.437 


-198.58 


-266.51 


-95.229 


150.39 


212.01 


111.12 


22.976 


- 


- 


- 


1 





1 


1 


-2.199 


-2.705 


-2.005 


-1.117 


0.402 


1.497 


2.562 


2.329 


1.827 


0.918 


0.424 


1 





- 1 


- 1 


4.011 
0.900 


2.215 
1.289 


3.875 
-0.409 


2.001 
0.691 


2.916 
0.750 


-0.045 
1.967 


0.723 
1.710 


-1.118 
1.584 


1.620 
0.746 


0.861 
0.307 


2.834 





1 


1 


X 


-2.798 
- 


0.641 
- 


-6.514 

- 


5.436 
- 


-11.528 


8.090 


-11.12 


5.779 

- 


-5.968 
- 


1.872 
- 


-1.377 

- 





1 


-1 


-1 


3.538 
0.448 


1.523 
0.479 


2.114 


0.383 


1.654 


-0.325 


1.1333 


0.130 


1.499 


0.773 


1.345 


1 


1 


1 


1 


62.71 

- 


417.41 

- 


1437.3 

- 


3228.7 

- 


5218.1 


6285.7 


5717.3 


3879.2 

- 


1894.7 

- 


607.32 

- 


101.69 

- 


1 


1 


-1 


-1 


88.678 


427.56 


1298 


2637.8 


4045.7 


4716.2 


4307.9 


2996.2 


1573.5 


555.82 


118.2 


Strategy 


a: 










Two Pietarinen expansions 










Ci 


Ca 


Bi 


B2 


Co 
do 


Cl 

di 


C2 
d2 


C3 

da 


C4 

dn 


C5 

d5 


C6 

de 


c? 
dr 


C8 
ds 


cg 
d(j 


610 

dio 


1 





1 


1 


3.849 


5.812 


-3.141 


-10.608 


-6.01 


-0.704 








- 




3.849 


-4.385 


4.217 


-3.331 


1.494 


-0.416 












1 





-1 


-1 


1.731 
1.731 


-2.290 
-1.212 


-0.002 
0.573 


-5.747 
-0.402 


0.693 
0.194 


-2.229 
-0.070 


1.239 
0.018 


- 


- 


- 


- 





1 


1 


1 


18.667 


117.70 


189.14 


171.86 


89.342 


22.696 




- 


- 


- 


- 


18.667 


-8.027 


8.294 


-5.408 


2.066 


-0.441 




- 


- 


- 


- 





1 


-1 


-1 


-32.588 
-32.588 


-217.23 
-0.676 


-220.56 
0.347 


30.972 
-0.195 


174.4 
0.109 


-5.098 
-0.055 


-190.34 
0.025 


-147.91 
-0.008 


-37.90 
0.002 


- 
- 


- 
- 


1 


1 


1 


1 


-1.255 


-2.575 


1.651 


4.307 


3.093 


-1.463 


-2.319 


0.232 


1.907 


1.221 


- 


-1.255 


-0.431 


-0.754 


1.189 


-1.122 


0.772 


-0.403 


0.144 


-0.032 


- 


- 


1 


1 


-1 


-1 


-14.963 
-14.963 


3.089 
136.57 


4.009 
-230.17 


3.351 
313.67 


3.146 
-443.08 


1.1584 
519.08 


0.1584 
-476.42 


-0.686 
301.89 


-0.572 
-87.63 


-0.430 

- 


- 
- 


Strategy 


a: 










Three Pietarinen expansions 


















CO 


Cl 


C2 


C3 


C4 


C5 


C6 


C7 


C8 


eg 


610 


Ci 


Ca 


Bi 


B2 


do 


di 


d2 


d3 


d4 


d5 




d7 


ds 


dg 


dio 










eo 


ei 


£2 


63 


£4 


£5 




67 


es 


eg 


eio 










2.275 


19.465 


27.86 


35.107 


18.749 


12.32 












1 


1 


1 


1 


2.274 


-0.616 


0.044 


-0.023 


0.019 


-0.007 




















2.277 


-0.878 


0.063 


-0.029 
























10.107 


1.372 


0.849 


-0.515 
















1 


1 


-1 


-1 


10.107 


-70.87 


65.949 


-21.659 
























10.107 


-0.757 


0.077 


-0.017 
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Pietarinen expansion coefficients corresponding to solutions presented in Table ^ (Strategy b). 



Startcgy b: 










One Pietarinen expansion 










Ci 


C2 


Bi 


B2 


CO 


ci 


C2 




C4 


C5 


C6 


C7 


eg 


C9 


ClO 










Cll 


Cl2 


Cl3 


Cl4 


Cl5 


Cl6 


Cl7 


Cl8 


Cl9 


e20 


e2i 


1 











0.676 


-0.795 


0.947 


-1.948 


2.450 


—3.236 


4.100 


-4.084 


4.034 


-3.702 


3.640 


-2.938 


2.575 


-1.934 


1.520 


—0.864 


0.553 


—0.248 


0.121 











1 








-0.362 


0.536 


-0.410 


0.211 


-0.170 


0.132 


-0.042 


0.074 


-0.011 


0.021 


- 








1 


1 


-0.560 


-1.205 


-2.397 


-1.364 


-0.517 


-1.170 


-0.612 


-0.193 


-0.199 


- 


- 








-1 


-1 


0.564 


-0.020 


1.994 


1.060 


0.168 


2.039 


0.397 


0.598 


1.357 






1 





1 


1 


1.255 


1.799 


2.300 


1.412 


1.418 


2.447 


1.143 


1.178 


1.134 


0.484 


0.300 


0.245 






















1 





- 1 


- 1 


0.777 


0.986 


1.832 


1.284 


1.641 


0.661 


0.807 


1.056 


0.709 


1.583 


0.617 










1.029 


0.987 


0.782 


1.281 


1.237 


0.316 


1.570 


1.216 


-0.781 


-0.934 


- 





1 


1 


1 


1.388 


1.753 


5.598 


6.395 


11.384 


10.090 


13.046 


9.156 


9.852 


3.338 


2.805 





1 


-1 


-1 


1.565 

— 


-0.031 

— 


-0.235 

— 


1.319 

— 


1.635 


1.836 


1.316 


1.179 

— 


0.392 

— 


0.381 

— 


— 
— 


1 


1 


1 


1 


-0.088 


1.193 


-0.502 


2.507 


-0.769 


3.256 


-1.326 


5.019 


-0.963 


2.993 


-0.645 




U.Uo / 




















1 


1 


-1 


-1 


4.426 


-1.681 


3.943 


-3.801 


3.305 


0.263 


-0.650 


6.418 


-4.716 


6.584 


-4.286 










2.423 


-1.484 


— 


— 








— 


— 


— 


— 


Strategy b: 










Two Pietarinen expansions 










Ci 


Ca 


Bi 


B2 


CO 
do 


ci 
di 


C2 
d2 


C3 
ds 


C4 
d4 


C5 
ds 


C6 

dg 


C7 
d7 


C8 
dg 


eg 
dg 


ClO 
dio 


1 





1 


1 


0.512 


1.761 


1.072 


-0.373 


1.552 


1.780 












0.512 


-0.467 


0.104 


-0.375 


0.148 


-0.039 




— 


— 


— 


— 


1 





-1 


-1 


0.301 


-2.718 


3.135 


0.373 


-3.393 


1.769 


2.913 


— 


— 


— 


— 










0.301 


0.455 


-2.287 


-1.714 


4.815 


-3.150 


0.837 


- 


- 


- 


- 





1 


1 


1 


0.793 


0.662 


0.629 


1.267 


1.042 


0.081 




- 


- 


- 


- 


0.793 


0.572 


-0.120 


0.035 


-0.011 


0.003 




- 


- 


- 


- 





1 


-1 


-1 


-1.649 


-12.558 


2.551 


20.807 


-4.393 


-17.838 


-14.688 


— 


— 


— 


— 










-1.649 


29.437 


8.936 


12.801 


12.373 


-25.197 


7.424 


— 


— 


— 


— 


1 


1 


1 


1 


6.410 


33.756 


55.745 


59.202 


44.944 


21.583 


6.587 


— 


— 


— 


— 


6.410 


1.223 


-0.903 


0.437 


-0.194 


0.060 


-0.171 


— 


— 


— 


— 


1 


1 


-1 


-1 


1.960 


-2.298 


2.642 


-2.579 


2.374 


-1.659 


1.027 


-0.424 


0.142 


— 


— 










1.960 


-0.817 


0.279 


-0.581 


0.392 


0.115 


-0.197 


0.068 


— 


— 


— 


Strategy b: 










Three Pietarinen expansions 


















CO 


ci 


C2 


C3 


C4 


C5 


C6 


C7 


eg 


eg 


ClO 


Ci 


C2 


Si 


B2 


do 


di 


d2 


ds 


d4 


ds 


de 


d7 


dg 


dg 


dio 










eo 


ei 


62 


63 


64 


es 


ee 


£7 


eg 


eg 


eio 










-1.608 


-5.203 


5.819 


1.520 


-3.035 


4.982 


5.711 


-0.305 








1 


1 


1 


1 


-1.608 


0.515 


0.252 


-1.565 


-2.579 


-2.229 


-1.062 


-0.267 
















-1.608 


1.518 


1.380 


-0.229 


-0.321 


0.411 


0.598 


0.282 
















1.842 


19.214 


40.571 


49.783 


42.630 


32.165 


6.188 


3.091 








1 


1 


-1 


-1 


1.842 


11.827 


24.661 


-54.662 


49.917 


-30.450 


12.023 


-2.268 
















1.842 


-5.396 


17.925 


-20.529 


11.084 


-2.477 
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